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1.0  INTRODUCTION 

The  failure  modes  associated  with  fractures  in  fiber-reinforced  composites  differ  considerably  from 
those  of  homogeneous  isotropic  materials.  The  failure  modes  in  composites  are  typically  in  the  forms  of 
transverse  cracking,  delamination,  fiber  breaks,  matrix  yielding,  matrix  cracking,  fiber  pull-out  and 
fiber/matrix  debonding.  As  a  result,  it  is  a  very  difficult  task  to  predict  the  fracture  process  and  behavior  of 
composite  materials  exactly. 

In  the  past,  a  great  deal  of  effort  has  been  expended  to  investigate  the  fracture  behavior  of 
composites.  A  number  of  fracture  mechanics  theories  have  been  proposed.  These  theories  have  been 
reviewed  and  presented  extensively  in  References  1  and  2. 

In  this  report  a  microscopic  theory  which  was  originally  proposed  by  Mar  and  Lin®  has  been  modified 
into  a  new  theory.  This  new  theory  links  the  inherent  flaw  concept®'®,  which  postulates  that  a  l  mage 
zone  exists  at  the  tip  of  crack,  and  the  theory  of  crack  singularities  at  a  bi-material  interface®'^®.  This 
combined  theory  can  be  used  to  predict  the  notched  strength  of  organic  and  metal  matrix  composites  with 
either  a  crack  or  a  hole. 

The  following  sections  of  this  report  describe  the  methodology  used  in  the  analysis,  assumptions 
taken,  analytical  results  obtained  and  conclusions  made. 
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2.0  THEORY 

The  stress  distribution  at  the  crack  tip  in  a  thin  plate  for  a  homogeneous,  isotropic  elastic  solid  in 
terms  of  the  coordinates  shown  in  Figure  1  is  given  by  equation  (1) 


Ox  “  On 


{ ^ 

2r 


V2, 


0/„  .  e  .  30 

cos-  1  -  sin-  sin— 


ao 


0/„  .  0  .  30 

cos-  1  +sin-  sin— 


ao 

txy-ON  2^ 


'/2, 


.0  0  30 

Sin-  cos^  cos-^ 


where 


(1) 


On  =  gross  nominal  stress. 
For  an  orientation  directly  ahead  of  the  crack  (0  =  0) 


Ox  •  Oy  -  On 


/  3o 


and  xxy  -  0 


(2) 


Irwin’ ^  pointed  out  that  equation  (1 )  indicates  that  the  local  stresses  near  a  crack  depend  on  the 
product  of  the  nominal  stress  and  the  square  root  of  the  half-flaw  length.  He  called  this  relationship  the 
stress  intensity  factor  K,  where  for  a  sharp  elastic  crack  in  an  infinitely  wide  plate,  K  is  defined  as 


K  “  On  Vac 


(3) 


In  the  approach  using  linear-elastic  fracture  mechanics  (LEFM),  K  is  a  material  parameter  and  may  be 
determined  from  tests. 

For  a  finite-width  plate,  equation  (3)  is  modified  to 

K  “  Yon  VSo  (4) 


Where  Y  is  a  parameter  that  depends  on  the  plate  and  crack  geometry. 

To  develop  a  similar  concept  for  composite  materials,  the  assumptions  of  references  3,  9  were 
adopted  in  this  paper;  i.e.: 
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On 


Figure  1 .  Model  For  Equations  For  Stresses  At  A  Point  Near  A  Crack. 
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•  The  existence  of  an  inherent  flaw  (also  called  a  damage  zone)  at  the  edge  of 
a  hole  or  at  the  tip  of  a  crack. 

•  Fracture  of  a  filamentary  composite  initiates  at  a  crack  lying  in  the  matrix 
material  at  the  interface  of  the  matrix/filament. 

•  A  laminate  fails  whenever  the  principal  load-carrying  laminae  fails.  This 
implies  that  cracks  in  the  matrix  perpendicular  to  the  0°  filaments  are  the 
triggering  mechanism  for  the  final  failure. 

Based  on  the  above  assumptions,  the  following  theories  are  developed: 

Using  the  same  concept  of  stress  intensity  factor  as  is  formulated  above  for  isotropic  materials,  a 
material  parameter  similar  to  K  is  defined  for  composite  material  as: 


K  -  Yon  (ao)"" 


(5) 


where  m  is  the  order  of  singularity  of  a  crack  whose  tip  is  at  the  interface  of  two  different  materials  as 
shown  in  Figure  2.  Calculations  for  determining  m  are  presented  in  Reference  18. 

Note  that  K  has  a  different  dimension  from  K.  (K  has  a  dimension  of  “stress  times  length  to  the  1/2 
power”,  while  K  has  a  dimension  of  "stress  times  length  to  the  m  power”.) 

Although  some  composite  materials  (such  as  polymeric  matrix  composites)  fail  in  a  brittle  manner, 
a  damage  zone  does  develop  which  is  analogous  to  the  plastic  zone  for  ductile  materials.  Using  this 
concept  in  conjunction  with  equation  (5)  yields: 

K  «  Yon  (ao  +  Co)”'  (6) 


where  Cq  is  defined  as  an  inherent  flaw  size.  The  term  inherent  flaw  size  is  used  since  unnotched 
strength.  Oq.  of  a  composite  laminate  is  given  by  equation  (6)  for  the  case  of  vanishing  a^, 


K-ao  (Co)”’ 


(7) 


where  Yq  is  the  correction  factor  for  infinite  plate. 

It  should  be  noted  that  Cq  does  not  physically  refer  to  an  inherent  crack,  but  a  characteristic 
dimension  of  damage  zone  at  the  tip  of  a  notch  or  crack  prior  to  ultimate  failure. 

The  question  we  may  ask  now  is  whether  K  and  C^  are  material  constants.  Before  we  reach  a 
conclusion,  certain  equations  are  helpful  in  answering  this  question. 


4 


NADC-881 18-60 


Substituting  equation  (7)  into  equation  (6),  and  after  somejnanipulations,  we  obtain  the  following 
important  equations  that  will  be  used  to  determine  parameter  K.  Cq  and  notched  strength  of  composite 
laminates 


Co  = 


ao 


Yo 


- 1 


N  ) 


(8) 


K  -  ao*  YgOo 


(9) 


On  J_  Cq 
Oo  Y  3o  Co 


(10) 


or 


1 


K 


(11) 


where 
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In  the  following  sections,  K  will  be  called  the  equivalent  stress  intensity  factor  for  composite  materials. 


On 


On 


Figure  2.  Crack  Normal  To  The  Bi-Material  Interface  With  Inherent  Flaw.  Co- 
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3.0  DETERMINATION  OF  EQUIVALENT  STRESS  INTENSITY  FACTOR,  K 
AND  INHERENT  FLAW  SIZE 

Reference  8  provides  extensive  fracture  test  data  of  boron/aluminum  laminates  with  various 
proportions  of  O'  and  ±45=  plies.  Hence,  this  test  data  will  be  used  to  characterize  the  fracture  behavior  of 
boron/aluminum  composite  laminates. 

3.1  EQUIVALENT  STRESS  INTENSITY  FACTOR.  K 

From  Appendix  A,  the  order  of  stress  singularity  at  the  boron/aluminum  interface,  m  equals  .347. 
From  equation  (9),  the  equivalent  stress  intensity  factor  for  boron/aluminum  composite  laminate  with 
center  crack  can  be  written  as  follows: 


K  -  ai*^ao\ 


Yon' 

Oo 


-288 


-  1 


-.347 


(12) 


Equation  (12)  is  used  to  characterize  the  critical  equivalent  stress  intensity  factor  of  boron/aluminum 
laminates  with  various  layups. 


By  using  the  fracture  test  results  from  reference  8,  equivalent  stress  intensity  factors  were  calculated 
from  equation  (12)  and  are  tabulated  on  Tables  1  through  4  for  boron/aluminum  composite  with  laminate 
constructions  [Olgi.  [02/±45]s,  (±45/021s  and  [0/±45]s  respectively.  Note  the  test  results  shown  in  Tables  1 
to  4  are  average  test  results  of  reference  8.  As  shown  in  Tables  1  to  4,  K  values  seem  to  be  a  material 
property  and  vary  with_different  laminate  orientations.  K  values  are  also  plotted  on  Figures  3  to  6.  Kavg  'S 
the  average  value  of  K  from  all  the  crack  sizes.  As  shown  in  the  figures.  K^vg  can  be  approximately 
treated  as  a  material  constant.  It  has  to  be  pointed  out  here  that  Kavg  was  obtained  by  averaging  three 
different  widths  of  plate.  For  w  =  101 .6mm  K  values  are  almost  the  same  for  different  2ao/w  ratios. 

The  detailed  calculations  of  K  are  shown  in  Appendix  B.1 . 
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Table  1 .  Equivalent  Stress  Intensity  Factor  For  [0]6T  B/AI  Composite. 


w 

(mm) 

Bo 

(mm) 

5 

■ 

K 

(MPa(mm) 

m 

Error 

% 

19.1 

.25 

.025 

1.001 

.818 

1125.0 

.878 

7.4 

19.1 

.65 

.068 

1.003 

.6845 

1136.0 

.760 

11.0 

50.8 

1.25 

.05 

1.001 

.592 

1166 

.663 

12.0 

50.8 

2.55 

.1 

1.006 

.555 

1455 

.546 

-1.6 

50.8 

.3 

1.06 

.4003 

1479 

.3705 

-7.5 

50.8 

12.7 

.5 

1.189 

.3098 

1518 

.279 

-10.0 

101.6 

2.55 

.05 

1 .001 

.580 

1455 

.549 

-5.3 

101.6 

5.1 

.1 

1.006 

.468 

1435 

.443 

-5.3 

101.6 

15.25 

.3 

1.06 

.3024 

1399 

.295 

-2.6 

101.6 

25.4 

.5 

1.189 

.232 

1430 

.221 

-4.7 

=  1360  MPa  (mm) 
oo=  1672  MPa 


Y- Vsec||^j 


8 


NADC-881 18-60 


Table  2.  Equivalent  Stress  Intensity  Factor  For  (O2/  ±45)s  B/AI  Composite. 


w 

(mm) 

3o 

(mm) 

■ 

K 

(MPa(mm)  ^^^) 

Q 

Error 

% 

19.1 

.25 

.025 

1.001 

.878 

652.6 

.891 

1.5 

19.1 

.65 

.068 

1.003 

.782 

681.9 

.782 

0. 

50.8 

1.25 

.05 

1.001 

.721 

740. 

.686 

-4.8 

50.8 

2.55 

.1 

1.006 

.571 

683.1 

.570 

0. 

50.8 

7.6 

.3 

1.06 

.395 

697.8 

.389 

-1.5 

50.8 

12.7 

.5 

1.189 

.274 

639.0 

.293 

IBSili 

101.6 

2.55 

.05 

1.001 

.579 

695.1 

.572 

-1.2 

101.6 

5.1 

.1 

1.006 

.464 

680.6 

.464 

0. 

101.6 

15.25 

.3 

1.06 

.321 

711.5 

.310 

-3.5 

101.6 

25.4 

.5 

1.189 

.229 

674.6 

.232 

1.6 

=  685.6  MPa  (mm) 
00  =  800.1  MPa 


w 
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Table  3.  Equivalent  Stress  Intensity  Factor  For  [±45/02]s  B/AI  Composite. 


w 

ao 

S 

Y 

K 

Error 

(mm) 

(mm) 

(MPa(mm)®^^) 

K|n 

% 

19.1 

.25 

.025 

1.001 

— 

_ 

19.1 

.65 

.068 

1.003 

_ 

_ 

50.8 

1.25 

.05 

1.001 

.619 

674.7 

.633 

2.3 

50.8 

2.55 

.1 

1.006 

.531 

716.6 

.517 

-2.6 

50.8 

7.6 

.3 

1.06 

.361 

720.8 

.349 

-3.5 

50.8 

12.7 

.5 

1.189 

.245 

647.7 

.262 

7.0 

101.6 

2.55 

.05 

1.001 

.528 

707.6 

.520 

-1.6 

101.6 

5.1 

.1 

1.006 

.423 

703.8 

.417 

-1.3 

101.6 

15.25 

.3 

1.06 

.281 

706.2 

.276 

-1.6 

101.6 

25.4 

.5 

1.189 

.201 

673.0 

.207 

3.1 

=  693.8  MPa(mm) 
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Table  4.  Equivalent  Stress  Intensity  Factor  For  [0/±45]s  B/AI  Composite. 


w 

(mm) 

a© 

(mm) 

K 

(MPa(mm) 

Q 

ERROR 

% 

19.1 

.25 

.025 

1.001 

_ 

— 

19.1 

.65 

.068 

1.003 

.828 

564.6 

.865 

4.4 

50.8 

1.25 

.05 

1.001 

.796 

645.8 

.789 

-1.0 

50.8 

2.55 

.1 

1.006 

.730 

710.1 

.688 

-5.7 

50.8 

7.3 

.3 

1.06 

.458 

597.8 

.481 

5.2 

50.8 

12.7 

.5 

1.189 

.329 

563.1 

.367 

11.5 

101.6 

2.55 

.05 

1.001 

.699 

656.5 

.692 

-1.0 

101.6 

5.1 

.1 

1.006 

.620 

702.0 

.569 

-8.2 

101.6 

15.25 

.3 

1.06 

.413 

678.0 

.388 

-6.0 

101.6 

25.4 

.5 

1.189 

.271 

583.9 

.293 

8.2 

=  633.5  MPa{mm) 
00  =  581.4  MPa 

Y- Vsec||^J 
-  _  2a 
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The  above  tests  are  for  center  crack  specimens.  For  other  crack  types  and  locations®,  the  calculated 
equivalent  stress  intensity  factors  are  shown  on  Table  5,  condensed  from  Appendix  B.1 .2.  It  can  be  seen 
that  K  for  these  unidirectional  boron/aluminum  composites  is  constant  for  different  crack  conditions.  The 
reason  for  the  slightly  different  K  as  compared  to  Table  1  is  due  to  a  different  value  of  ultimate  tensile 
strength. 


Table  5.  Critical  Equivalent  Stress  Intensity  Factor  Of  Unidirectional  B/AI®  Composite. 


SPECIMEN* 

2ao  2R 

—  or  — 

Kavg  MPa  (mm) 

w  w 

CH 

.25 

1246.00 

CS 

.40 

1227.00 

DEN 

.30 

1237.00 

*CH  -  CENTER  HOLE  SPECIMEN 
CS  -  CENTER  SPLIT  SPECIMEN 
OEN  -  DOUBLE  EDGE  NOTCH  SPECIMEN 

3.2  INHERENT  FLAW  SIZE,  Cq 

Two  methods  were  used  to  calculate  the  inherent  flaw  sizes  for  a  composite  laminate  with  center 
crack, 

Least  Square  Fit 

Equation  (8),  in  which  Cq  is  a  proportional  constant,  can  be  rearranged  to  yield 


-2  88  1 
-1 

|i  O,  ; 

! 

) 

By  using  the  fracture  test  data  as  shown  in  Tables  1  to  4,  and  the  least  square  fit,  Cq  for  various 
laminate  constructions  are  determined  as  shown  in  Table  6. 

Average  Equivalent  Stress  Intensity  Method 

From  Equation  (7)  we  have 

KAVG  =  ao(Cor  (14) 

where  m  =  .347  for  boron-aluminum. 

The  inherent  flaw  size  can  be  derived  from  equation  (14)  as  follows: 
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^Kavg^ 

Oo 


(15) 


In  the  case  of  boron/aluminum  composite,  equation  (15)  becomes 


Co- 


,2.88 

(Kavg\ 


Oo  J 


(16) 


The  inherent  flaw  sizes  for  various  laminate  orientations  were  calculated  using  this  method  and  were 
also  tabulated  on  Table  6. 


Table  6.  Inherent  Flaw  Size,  Co  (mm). 


Least  Square  Fit 

Kavg  Method 

l'0]6T 

.633 

.552 

(02/:t45)s 

.593 

.641 

[±45/0a]s 

.405 

.457 

[0/±45]s 

.976 

1.28 

We  pointed  out  earlier  that  Co  does  not  physically  refer  to  an  inherent  crack,  but  rather  to  a 
characteristic  dimension  of  damage  zone  at  a  crack  tip,  prior  to  fracture  failure.  Comparing  the  dimension 
of  Co  with  respect  to  the  crack  size  of  B/AI  specimens  as  shown  in  Table  1  to  4,  it  is  clear  that  Co  cannot 
be  neglected  in  the  calculations  of  equivalent  stress  intensity  factor.  The  significance  of  Co  will  be  further 
discussed  in  Section  4.0. 

Equation  (7)  is  used  to  calculate  critical  equivalent  stress  intensity  factor  based  on  the  Co  determined 
from  least  square  fit  method.  These  values  are  also  plotted  on^Figures  3  to  6  as  Klsf-  It  can  be  seen  from 
the  plots  that  there  is  not  much  difference  between  Rlsf  and  Kavg  excepHn  th^plot  for  the  [0/s45]s 
laminates,  where  though  there  is  a  greater  difference  between  Kavg  and  Klsf,  I^vg  provides  a  better 
result.  For  this  reason  Kavg  will  be  adopted  in  this  report  and  will  be  denoted  as  Kq. 
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Figure  3.  Equivalent  Stress  Intensity  Factor  For  [0]6T  B/AI  Composite. 


K 

MPa  (mm) 


w 

Figure  4.  Equivalent  Stress  Intensity  Factor  For  [02/i:45]s  B/AI  Composite. 
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Figure  5.  Equivalent  Stress  Intensity  Factor  For  [±45/02ls  B/AI  Composite. 


,05  .1  .2  Pa  .3  .4  .5 


W 

Figure  6.  Equivalent  Stress  Intensity  Factor  For  [0/i:45]s  B/AI  Composite. 
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4.0  APPLICATIONS 


Once  the  critical  equivalent  stress  intensity  factor.  Kq  is  known,  the  fracture  strength  of  the  composite 
laminates  can  be  obtained  from  equations  (10)  and  (1 1). 

From  equation  (11 ),  we  have  a  fracture  strength  prediction  formula  as  follows: 


Y 


1  + 


Kr 


J 


(17) 


,  the  lesser  the  notch  sensitivity  and 


/Kq' 

From  Equation  (1 7),  it  can  be  seen  that  the  larger  the  term  — 

SO  from  equation  (7),  we  can  conclude  that  the  larger  the  inherent  flaw  size  (i.e.,  the  damage  zone  size), 
Co,  the  lesser  the  notch  sensitivity. 


In  the  following  subsections,  the  theory  developed  here  will  be  used  to  predict  the  fracture  strength  of 
various  composite  laminates. 

4.1  BORON/ALUMINUM  (B/AI)  COMPOSITE 

4.1.1.  Notched  Strenqth  Prediction 

For  B/AI  composite  laminates,  m  =  .347  and  for  a  center  crack  specimen,  equation  (17)  becomes 


On  1 
oo  Y 


_  ,-2.88  .  347 

ao] 

Oo 


(18) 


where  for  a  composite  laminate  with  a  center  crack,  Y  is  assumed  to  be  the  same  as  that  for  an  isotropic 
material®. 


For  convenience,  data  from  Tables  1  to  4  and  6  are  summarized  on  Table  7  to  be  used  for  the 
following  analysis.  Substituting  oo  and  Kq  from  Table  7  for  different  ply  conditions  into  equation  (18) 

values  for  —  can  be  obtained  for  various  crack  sizes  and  are  tabulated  on  Column  7  of  Tables  1  to  4.  The 

Oo 

analytical  results  are  also  plotted  on  Figures  7  to  10. 

As  can  be  seen,  the  prediction  represents  the  experimental  results  reasonably  well. 

The  Kq  values  obtained  from  center  cracked  specimens  are  also  applied  to  B/AI  specimens  with 
center  holes.  Comparisons  of  the  analytical  results  and  test  results^'*  are  shown  in  Figures  1 1  to  13. 

Figure  1 1  shows  excellent  correlation  between  test  and  analytical  results  for  [0]6T  laminates  with  center 
holes,  while  in  Figure  12  and  13,  the  meWimum  percentage  error  for  analytical  results  is  around  13%.  The 
detailed  calculations  are  shown  in  Appendiz  B.1  3.  This  confirms  the  findings  of  reference  18  and  19 
that  the  length  of  discontinuity  and  not  the  shape  appeared  to  control  the  fracture  strength  of 
composite  laminates. 
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Table  7.  Fracture  Parameters  For  Various  Laminate  Configurations  Of  B/AI. 


Ply 

Configuration 

Oo 

(MPa) 

Kq 

MPa  (mm) 

Oo 

(mm)3"7 

Co 

(mm) 

[0]6T 

1672 

1360 

.81 

.552 

[02/±45Is 

800.1 

685.6 

.867 

.641 

[i:45/02]s 

910.5 

693.8 

.762 

.457 

[0/445]s 

581.4 

^  -  -  — - 

633.5 

1.09 

1.28 
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Notched  Strength  Results  For  B/AI  [02/i45]s  Composites. 
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e  9.  Comparison  Of  Predicted  And  Experimental  Notched  Strength  Results  For  B/AI  (±45/02]s  Composites. 
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Figure  10.  Comparison  Of  Predicted  And  Experimental  Notched  Strength  Results  For  B/AI  (0/±45]s  Composites. 
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Figure  1 1 .  Notched  Strength  Prediction  For  B//AI  (OJer 
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Figure  12.  Notched  Strength  Prediction  for  B/AI  (02/±45]s 
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Figure  13.  Notched  Strength  Prediction  For  B/AI  [0/±45Js 
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4.1.2.  Notch  Sensitivity  of  Boron/Aluminum  (B/AI)  Composite  Laminate 


The  test  data  in  Tables  1  to  4  are  plotted  on  Figures  14  to  16  for  composite  laminate  of  various 

Oo 

widths  to  check  the  notch  sensitivity  of  B/AI  composites.  It  is  obvious  that  [±45/02]s  is  the  most 

23o  -  ■ 

notch-sensitive,  while  [0/±45]s  is  the  least  sensitive  to  the  notch  size.  For  —  s  .1 ,  [OJst  is  more  sensitive 


to  the  notch  size  than  [Oa/^^SJs.  For  >  .1 ,  the  ^  vs  curve  for  [0]6T  and  [02/±45]s  laminates  are 
almost  identical.  The  trend  mentioned  above  can  be  detected  by  using  a  parameter^®  defined  by  the  ratio 
~  as  shown  in  Table  7.  The  ranking  of  various  laminate  configurations  for  the  —  index  is  [0/±45]s, 

Oo  Oo 

r7_ 

[02/2:45]s,  [0]6t  and  [s:45/02]s  in  descending  order.  The  relative  values  of  —  show  a  similar  trend  to  that  of 

Oo 

damage  zone  size,  Co¬ 


lt  can  be  concluded  that  the  larger  the  ratio 
sensitivity  there  is  to  the  notch  size. 


kq 

Oo 


(or  the  damage  zone  size  Co)  the  less  fracture 


4  2  GRAPHITE/EPOXY  (Gr/Ep)  COMPOSITE 

Equation  (9)  is  used  Jo  characterize  the  Rq  for  Gr/Ep  composites.®  The  detailed  calculations  are 
shown  in  Appendix  B.2.  Kq  for  various  ply  orientations  is  plotted  as  shown  in  Figures  17  to  19.  It  can  be 
seen  that  Kq  for  Gr/Ep  composites  can  be  treated  as  a  material  constant.  Table  8  summarizes  the 
characterization  results  of  Gr/Ep  composite  laminates.  For  Gr/Ep.  m  =  .297. 
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4.2.1  Notched  Strength  Prediction  and  Notch  Sensitivity 
For  Gr.  Ep  composite  with  center  crack,  equation  (17)  becomes 


On 

Oo 


1 

Y 


1  + 


.  —  .  -3.367  _  297 

(ooj 


(19) 


Substituting  Ko  and  oo  from  Table  8,  into  equation  (19),  the  fracture  strengths  of  graphite/epoxy  for 
various  ply  orientations  can  be  obtained  and  are  plotted  on  Figure  20.  The  detailed  calculations  are 
shown  in  Appendix  B.2.  As  can  be  seen  from  the  figure,  the  correlation  between  analytical  and 
experimental  results  is  very  good. 


It  can  also  be  seen  from  Figure  20  that  the  ratio  —  (or  the  inherent  flaw  size)  in  Table  8  can  be  used 

Oo 

as  a  notch  sensitivity  indicator.  The  [0/90/±45]s  laminate  is  less  notch  sensitive  than  the  (0/±45]s  and 
[0/±45]2s  laminates  and  accordingly  it  has  a  larger  ratio  of  --  (or  Co)  than  the  other  two  laminates. 

Oo 
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Table  8.  Fracture  Parameters  For  Various  Laminate  Configurations  Of  Gr/Ep. 


Ply 

Configuration 

Oo 

(MPa) 

Kq 

MPa  (mm)'^^ 

Oo 

(mm)297 

Co 

(mm) 

[0/±45]2s 

541.0 

408.6 

.755 

.389 

(0/*45]s 

541.0 

393.5 

.727 

.342 

[0/90/±45]s 

454.0 

437.7 

.964 

.884 
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4.3  GLASS/EPOXY  (Gl/Ep)  COMPOSITE 

Equation  (9)  is  also  used  to  characterize  the  Rotor  glass/epoxy  composites. ’6  The  detailed 

23 

calculations  are  shown  in  Appendix  B.3.  Rq  for  [0'±45/90l2s  ply  orientation  is  plotted  against  —  in 
Figure  21 .  For  Gl/Ep  composite  m  =  .289.  As  seen  in  the  figure,  Rq  can  be  approximately  treated  as  a 
material  constant. 


For  GICp  composites  with  center  crack,  equation  (17)  becomes 


with 


On 

Oo 


1 

Y 


1  + 


,-3.46 

IS]  ‘I 


Kq  =  274  MPa  (mm) 
Oo  =  320  MPa 


Equation  (20)  becomes 

On  1  a  |-.289 

---1  +  1.7108  aol 

Oo  Y  ' 

For  the  laminate  with  a  center  crack®, 

where  ao  =  half  crack  length. 


(20) 


(21) 


(22) 
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For  the  laminate  with  a  center  hole.  Y  is  assumed  in  the  following  form 


Y  = 


(23) 


where  Ro  =  radius  of  the  hole 

Equation  (21 )  was  plotted  for  Gl/Ep  [0/±4S/90]2s  laminates  with  a  center  crack  and  a  hole  as  shown  in 
Figure  22. 

It  can  be  seen  that  the  correlations  between  analytical  results  and  test  results  are  good.  Also,  note 
that  a  center  hole  decreases  the  fracture  strength  of  a  composite  laminate  slightly  more  than  a  center 
crack  does. 


4.4  COMPARISON  OF  ANALYTICAL  RESULTS  BETWEEN  ANISOTROPIC  AND  MICROSCOPIC 

MODELS 

Composite  materials  made  by  combining  two  materials  with  different  elastic  modulii  are  by  nature 
anisotropic  in  the  gross  sense.  The  anisotropic  model  for  composite  materials  is  to  assume  that  the 
composite  is  a  homogeneous,  anisotropic  solid.  For  an  anisotropic  fracture  m  =  .5.^^  Applying  the 
inherent  flaw  concept  for  an  anisotropic  model,  for  center  crack  specimen,  we  have 

Kq  -  YoN(ao  CS)’'’  (24) 

Kq  -  Oo(C5)  (25) 

Note  that  Kq  has  dimensions  different  from  those  of  Kq  and  Co  is  the  inherent  flaw  size 
corresponding  to  an  anisotropic  model. 
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Figure  22.  Fracture  Strength  Prediction  For  Gl/Ep  [0/±45/90]2s. 
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Examining  equations  (24)  and  (25),  we  can  derive  the  following  useful  equations  for  an  anisotropic 
model 


CS- 


a© 


On 

Oo 


1 

Y 


yBo  +  Coj 


(26) 


(27) 

(28) 


Equation  (26)  can  be  used  to  obtain  CS  through  the  least  square  fit  method.  It  has  been  determined 

for  the  anisotropic  model  that  Co  found  by  the  least  square  fit  method  predicts  the  test  data  better  than  Co 
found  by  the  average  equivalent  stress  intensity  method.  Equation  (27)  can  be  used  to  predict  the  fracture 
strength  of  composite  laminates. 


In  this  report,  only  B/AI  composites  will  be  used  to  demonstrate  the  superiority  of  the  microscopic 
model  over  the  anisotropic  model. 


Tables  9  to  1 2  show  the  comparison  of  analytical  results  between  the  microscopic  and  anisotropic 
models.  These  results  are  also  plotted  on  Figures  23  to  26.  It  is  clear  that  the  microscopic  model  predicts 
better  results  than  the  anisotropic  model.  Note  that  the  results  of  the  microscopic  model  are  copied  from 
Tables  1  to  4,  while  the  results  of  the  anisotropic  model  are  shown  in  Appendix  C. 
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Figure  23.  Comparison  Of  Analytical  Results  Between  Anisotropic  And  New  Approximate  Models  B/AI  (O)eiT- 


Anisotropic  Models 
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Figure  24.  Comparison  Of  Analytical  Results  Between  Anisotropic  And  New  Approximate  Models  B/AI  (02/+45)s. 


Comparison  Of  Analytical  Results  Between  Anisotropic  And  New  Approximate 


TEST 
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Figure  26.  Comparison  Of  Analytical  Results  Between  Anisotropic  And  New  Approximate  Models  B/AI  (0/+45), 


Table  9.  Comparison  Of  Analytical  Results  Between  Anisotropic  And  New  Approximate  Models  B/AI  (Olej. 
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Table  10.  Comparison  Of  Analytical  Results  Between  Anisotropic  And  New  Approximate  Models  B/AI  [Oj/tAS), 
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Table  11 .  Comparison  Of  Analytical  Results  Between  Anisotropic  And  New  Approximate  Models  B/AI  [+45/O2], 
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Table  12.  Comparison  Of  Analytical  Results  Between  Anisotropic  And  New  Approximate  Models  B/AI  (0/±45), 
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5.0  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  CONCLUSIONS 

•  The  methodology  developed  here  can  be  used  to  characterize  the  fracture  toughness  of  the 
composite  laminates  and  can  be  used  as  a  design  tool  to  predict  the  fracture  strength  of  various 
composite  laminates. 

•  The  parameter  Kq  which  was  called  critical  equivalent  stress  intensity  factor  is  defined,  and  can 
be  treated  as  a  material  constant  for  composite  laminate. 

•  The  new  approximate  model  provides  better  results  than  those  of  the  anisotropic  model. 

•  The  larger  the  ratio  —  (or  the  inherent  flaw  size.  Cq).  the  higher  the  damge  tolerance. 

Oo 

5.2  RECOMMENDATIONS 

•  Further  verification  of  microscopic  theory  with  test  results  of  various  composites  materials  is 
needed. 

•  Apply  the  theory  developed  here  to  predict  the  fracture  strength  of  composite  laminates  with 
various  crack  angles. 

•  Develop  a  methodology  to  predict  the  inherent  flaw  size  at  the  crack  tip  before  fracture. 
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APPENDIX  A 

CALCULATION  OF  THE  ORDER  OF  STRESS 
SINGULARITY  AT  A  BI-MATERIAL  INTERFACE 

In  the  case  of  a  crack,  normal  to  the  bi-material  interface  as  shown  in  figure  A-1 ,  the  characteristic 
equation  to  determine  the  order  of  stress  singularity  is  given  as  follows.’®; 

+  43cp)  +  23c^-2ap  +  23t-p  +  1+(-2=c2+ 2ap-2a  +2p)coST:t  =  0  (A-1) 

where 


P 


i^-1 

,  _  F2 _ 

1 

FI  (h2  1) 


F2(ril  +  1) 

T11  -  3  -  4vi  ,  ti2  -  3-  4v2  for  plane  strain 
ni  -  (3  -  V2)/(1  +  vi) ,  ti2  -  (3  -  V2)/(1  +  V2)  for  plane  stress 
FI ,  F2  =  Shear  modulus  of  medium  1  and  2  respectively 
VI,  V2  =  Poison’s  ratio  for  medium  1  and  2  respectively 
The  stresses  near  the  crack  tip  (for  0  =  0)  can  be  written  as 

.T-  1 


Ox  -  r 


Oy  -  r 


T-  1 


Xxy-  r 


T-  1 


The  order  of  stress  singularity  is  defined  as 


m  =  1  -T 


(A-2) 


(A-3) 


(A-4) 


By  treating  the  matrix  as  medium  1 ,  and  fiber  as  medium  2,  equations  (A-1)  to  (A-4)  are  used  to 
calculate  the  order  of  singularities  for  various  composite  materials  as  follows; 

BORON/ALUMINUM 


The  properties  of  boron  fiber  and  aluminum  matrix  are  as  follows; 
Aluminum;  m  =  3.76  msi,  vi  =  .33 

Boron:  H2  =  26.77  msi,  V2  =  .13 
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Internee 


Figure  A-1 .  Crack  Normal  To  The  Bi-Material  Interface. 


Substituting  these  properties  into  equations  (A-1)  and  (A-2)  for  plane  stress  case,  we  have  the  following 
characteristic  equation; 

.5191  -  .6448  cos  T  n  -  .5204  =  0  (A-5) 

Solving  equation  (A-5),  we  have 

T  =  .653 

From  equation  (A-4) 

m  =  .347 

GRAPHITeEPOXY 

For  Thornel  graphite  fiber  properties 

H2  =  4  msi  V2  =  .2 

For  Epoxy  matrix  properties 

pi  =  ,19  msi  VI =.35 

Applying  the  same  procedures  as  for  the  Boron/Aluminum  Composite,  we  find 

m=.297  for  Thornel  Qraphite/Epoxy 
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GLASS/EPOXY 

For  E-type  glass  fiber, 


For  Epoxy  matrix, 


(12  =  4.4  msi  V2  =  .2 


m  =  .17msi  vi  =  .35 

Applying  the  same  procedure  as  for  the  Boron/Muminum  composite  we  find: 

m  =  .289  for  E-type  Glass/Epoxy 
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APPENDIX  B 

CALCULATION  OF  Kq  AND  —  FOR  VARIOUS  COMPOSITE  LAMINATES 

Oo 

The  following  equations  were  used  to  calculate  K  and  ~ 


K 


^  fri  „ 

ao  Oo 


Oo 


1 

Y 


f  1 

1  + 1 

Ool 

3o 

(B-1) 


{B-2) 


B.1  BORON/ALUMINUM  COMPOSITE  (M  =.347) 
B.1 .1  Determination  of  Kq  fT ables  1  to  4) 


Column  5  of  Tables  1  to  4  listed  the  test  results  of  Substituting  these  quantities  into  equation 

Oo 

—  2ao  — 

(B-1),  we  obtain  K  for  various  S  (i.e.  —  )  values  as  shown  in  Column  6  of  Tables  1  to  4.  Kq  is  then 

defined  as 


n 

1  V*  75 


Kq  -  KaVG  -  ^  ^ 


i-  1 


(B-3) 


Where  n  =  total  number  of  test  data,  and  where  Kavg  rather  than  Klsf  is  chosen  to  be  Kq.  the  critical 
stress  intensity  factor  as  demonstrated  in  Figures  3  through  6  in  Section  3.2. 

B.1 .2  Determination  Of  Kq  From  Test  Data  Of  Reference  5 


Column  4  of  Tables  B-1  to  B-3  shows  the  test  results  for_yarious  ratios  of  Using  the  same 
procedures  described  in  the  previous  section  we  can  obtain  Kq  for  composite  laminates  with  a  center 
hole,  double  edge  notch,  and  center  slit  respectively  as  shown  in  Column  5  of  Table  B-1  to  Table  B-3. 

Table  B-1 .  Unidirectional  B/Ai  Composite  Laminates  With  Center  Hole. 


oo  =  1470.79  MPa 


Ro 

(mm) 

S 

Yh 

— 

R 

MPa  (mm)-®^^ 

On 

Oo 

1 .5875 

.25 

1.04 

.542 

1048.31 

.601 

2.38125 

.25 

1.04 

.521 

1149.86 

.538 

3.175 

.25 

1.04 

.513 

1247.19 

.495 

6.35 

.25 

1.04 

.443 

1339.3 

.400 

Kq  =  Kavg  =1196  MPa  (mm)  3^^ 
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Table  B-2.  Unidirectional  B/AI  Composite  Laminates  With  Double  Edge  Notch. 


ao 

(mm) 

fc 

s. 

K 

MPa(mm)^^ 

On 

Oo 

1.905 

.3 

1.01 

.565 

1270.02 

.546 

2.8575 

.3 

1.01 

.453 

1128.47 

.489 

3.81 

.3 

1.01 

.471 

1302.47 

.443 

7.62 

.3 

1.01 

.338 

1156.97 

.354 

Kq  =  Kavg  =  1214.5  MPa  (mm) 

Yc  =  (1 .98  +  .36  E  -2.12 +  3.42  £^)/1 .77 
Table  B-3.  Unidirectional  B/AI  Composite  Laminates  With  Center  Slit. 


ao  1 

Yh 

(-1 

K 

On 

(mm) 

l°°iTEST 

MPa  (mm) 

Oo 

3.81 

.4 

1.11 

.431 

1160.85 

.454 

7.62 

.4 

1.11 

.382 

1292.61 

.362 

Kq  =  Kavg  =  1227  MPa  (mm) 


Yc-Vsec|^j 

Table  B-4.  Fracture  Strength  Prediction  For  B/Al  [Ojer  With  Center  Hole. 


Ro 

(mm) 

I  1 

I  Yh  I 

l  j 

Rq  .347  1 

MPa(mm) 

i 

ft 

1  /  V 

1  1—^ 
(ooj 

1  Error 

% 

1 

1.59 

1.0625 

!1.0  1 
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1 
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il.01  i 
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:  -5.0 

6.35 

.125 

1.01 

1360  1 
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^  .412 

-1.7 

6.35 

.25 

!  1 .04  i 
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.4 

J _ 

i 

12.7 

.25 

1  i 

i1.04 
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1360  i 
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1 

'  .319 

1 

'  6.7 

—  (1  +  1.813Ro)-^^^ 

Oo  » 
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B.1 .3  Prediction  of  —  for  Boron/ Aluminum  Laminates  with  Center  Hole 
- oo - 

From  Tables  1  to  Table  4,  we  have  Kq  for  [0]6T,  t02/±45]s  and  [0/±45]s  as  1360  MPa  (mm)  685.6 
MPa  (mm)  and  633.5  MPa  (mm)  respectively.  Substituting  these  quantities  into  equation  (B-2)  for 
K,  we  can  predict  the  fracture  strength  of  composite  laminates  with  a  center  hole  as  shown  in  Tables  B-4 
to  B-6.  The  test  data  are  obtained  from  Reference  14. 


Table  B-5.  Fracture  Strength  Prediction  For  B/AI  {02/±45]s  With  Center  Hole. 


Ro 

(mm) 

t 

•"I 

Yh 

! 

Kq 

P). 

H 

Error 

% 

1.59 

'  .0625 

1.0 

! 
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.649 
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'  .125 

1.01 

685.6 
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.533 

1.5 

16.35 

i.125 

1.01 

i 
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_ 1 
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i  .432 

-9.2 

6.35 

,.25 

1.04 

685.6 

.450 

.420 

-6.7 

12.7 

1.25 

1.04 

685.6 

.375 

!  .335 

10.0 

On 

0  +  1.56RJ 

-.347 

Oo 

Yh 

Table  B-6.  Fracture  Strength  Prediction  For  B/AI  [0/±45]s. 

ao 

£ 

Yh 

! 

Kq 

(-1 

I  (On\ 

Error 

(mm) 

1  looj 

! 

% 

1.59 

.0625 

1.0 
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'  .756 

—4.6 
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1.125 

1.01 
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.645 

.642 

0. 

6.35 

1.01 

1 
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.585 
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-9.0 
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i.25 
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1 
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♦ 
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i  .419 

-10.0 

— -;i(1  +  .781Ro)"®^^ 

Oo  • 
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B.2  GRAPHITE/EPOXY  COMPOSITE  (m=.297) 

Substituting  m=.297  and  the  t^t  data  from  Reference  6  as  shown  in  Column  4  of  Tables  B-7  to  B-9, 
into  equation  (B.1),  we  can  obtain  K  as  listed  on  Column  5  of  the  Tables.  Then  use  equation  (B-3)  to 
obtain  Kq. 

After  obtaining  Kq,  substituting  these  quantities  into  equation  (B-2),  we  can  obtain  —  of  various  £  for 

Oo 

different  laminate  configurations  as  shown  in  Column  6  of  Tables  B-7  to  B-9. 

B.3  GLASS/EPOXY  COMPOSITES  (m=.289) 

Substituting  m=.289  ^d  the  test  data  from  Reference  16.  as  shown  in  column  4  of  Table  B-10,  into 
equation  (B-1)  we  ootain  K  as  listed  in  column  5.  Kq  can  be  obtained  from  equation  (B-3).  By  substituting 

Kq  into  equation  (B-2),  we  can  obtain  —  for  the  composite  laminate  [0/±45/90]s  with  both  a  center  crack 

Oo 

and  a  center  hole  as  shown  in  Column  6  of  Tables  B-10  and  B-1 1 . 


Table  B-7.  Fracture  Parameters  Of  Gr/Ep  Laminate  [0/±45]2s  With  Center  Crack. 


a© 

(mm) 

[ 

■ 

Test  ! 

R  ' 

[MPa  mm) 

» 

H 

2.54 

.1 

11.01 

.56 

i 

423 

r 

1 

( 

.543 

m 

5.08 

.2 

1.026 

.46 

i 

424 

.444 

-3.4 

7.62 

:  .3 

^  1.054 

.38 

i 

1 

402 

1 

i 

.386 

1.6 

10.16 

.4 

1.103 

.34 

1 

i 

409 

i 

i 

.34 

11 

12.7 

.5 

1.189 

.28 

386 

i 

.296 

m 

Rq  =  Kavg  =  408.6  MPa(mm) 


Oo  =  541  MPa 
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Table  B-8.  Fracture  Paranrreters  Of  Gr/Ep  Laminate  [0/*45]s  With  Center  Crack. 


ao 

(mm) 

'  £ 

1  ■ 

1 

■ 

m 

!  K 

1  (MPa  mm) 

/On) 

Error 

% 

1 

i 

mm 

i 

11.01  1 

.53 

397 

.526 

0. 

5.08 

i 

1.026  1 

.44 

404  1 

.429  ! 

-2.5 

7.62 

!  -3 

1 1.054  ! 

\  \ 

.39 

1 

413 

1  .373 

-4.5 

10.16 

'  .4 

!  1.103  ! 

.33 

396 

1  .328 

i  -.6 

12.7 

'  .5 

1.89  ■ 

.26  ' 

'  358 

.285 

9.7 

Ko=Kavg=393.5  MPa{mm) 


oo  =  541  MPa 


Table  B-9.  Fracture  Parameters  of  Gr/Ep  Laminate  [0/90/±45]s  With  Center  Crack. 


ao 

(mm) 

1 

i 

1 

Y  ' 

1 

i  i 

m 

K 

(MPa  mm) 

H 

1 

!  Error 

1 

2.54 

1 

i  .1 

j _ 

li.oi  i 

.69 

1 

464 

i 

.662 

i 

5.08 

i  .2 

i 

1.026 

.57  ! 

'  454  ' 

.552 

-3.1 

7.62 

.3 

1.054 

.47 

424 

.484 

2.9 

!  .4 

,  1.103  1 

.416 

1 

425 

.428 

i  2.8 

12.7 

1 

,  .5 

i  1.189  : 

.36 

422 

.373 

1 

'  3.7 

Rq  ='Kavg  =  437.7  MPa(mm) 
Oo  =  454  MPa 
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Table  B-10.  Fracture  Parameters  of  Gl/Ep  (0'±45/90]s  with  a  Center  Crack. 


ao 

(mm) 

;  1  " 

1  ; 

a. 

1  R 

-  (MPa  mm) 

1  f£^) 

lool 

j  ' 

Error 

% 

1.45 

1  /II 

11.01  1 

.707 

1  283.3 

i  .708 

0. 

3.86 

/.30 

1.054  : 

.578 

!  305 

.530 

-6.4 

7.37 

I  .29 

i1.05  i 

.440 

1  268.8 

.448 

1.0 

12.4  1 

'  .33 

1.06 

.338 

239.4 

.382 

13.0 

Kq=  Ravg  =  274  MPa(mm) 
Oo=320MPA 


I 

Table  B-11 .  Fracture  Strength  Prediction  of  Gl/Ep  [0/±45/90]s  with  a  Center  Hole. 


Ro 

(mm) 

i 

i 

I  Y„ 

Rq 

:  (MPa  mm)  , 

(5) 

!  Error 

i  % 

1 

1.27 

1  .1 

1.01 

274 

.696 

.709 

1.8 

_ 

3.81 

:  .3 

:  1 ,054 

1  274 

.511 

.530 

1  3.6 

7.62 

.3 

1.054 

274 

.44 

.442 

-5 

12.7 

.33 

1.06 

274 

.375 

.383 

m 
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ANALYTICAL  RESULTS  OF  ANISOTROPIC  MODEL 
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APPENDIX  C.  ANALYTICAL  RESULTS  OF  ANISOTROPIC  MODEL. 


INHERENT  FLAW  CONCEPT 
B/AL  (±45/02)s  exponent  =  0.500 
UNNOTCHED  STRENGTH  =  910.5  Co  =  1 .505519 


a 


2a 

w 


Y 


Ooj 


exp-1 


/On' 

TEST 


On 
Oo 

CALC 


%ERR 


1.25 

0.05 

1 .00579E+00 

1.57989E+00 

6.19000E-01 

7.34908E-01 

18.7 

2.55 

0.05 

1.00590E+00 

2.54504E+00 

5.28000E-01 

6.05709E-01 

14.7 

2.55 

0.10 

1.011 51 E+00 

2.46633E+00 

5.31000E-01 

6.02351  E-01 

13.4 

5.10 

0.10 

1.011 51 E+00 

4.46234E-I-00 

4.23000E-01 

4.71975E-01 

11.6 

7.60 

0.30 

1 .05343E+00 

5.91474E+00 

3.61000E-01 

3.85999E-01 

6.9 

12.70 

0.50 

1.18275E+00 

1.09092E+01 

2.45000E-01 

2.75246E-01 

12.3 

15.25 

0.30 

1 .05379E+00 

1.04046E+01 

2.81000E-01 

2.84453E-01 

1.2 

25.40 

0.50 

1.18275E+00 

1.66938E-»^01 

2.01000E-01 

2.00000E-01 

0.5 

INHERENT  FLAW  CONCEPT 
B/AL  (0/±45)s  EXPONENT  =  0.500 
UNNOTCHED  STRENGTH  =  581 .4  Co  =  2.836446 


2a 

w 


Y 


exp-1 


On 

Oo 


%ERR 


TEST  CALC 


0.65 

0.07 

1.00787E+00 

4.3591 4E-01 

8.28000E-01 

8.94933E-01 

8.1 

1.25 

0.05 

1 .00579E+00 

5.601 15E-01 

7.96000E-01 

8.28334E-01 

4.1 

2.55 

0.05 

1.00590E+00 

1.08791  E+00 

6.88000E-01 

7.21407E-01 

4.9 

2.55 

0.10 

1.01 151  E+00 

8.34061  E-01 

7.30000E-01 

7.17407E-01 

1.7 

5.10 

0.10 

1.01 151  E+00 

1.54259E+00 

6.20000E-01 

5.91022E-01 

4.7 

7.60 

0.30 

1.05343E+00 

3.29595E+00 

4.58000E-01 

4.94887E-01 

8.1 

12.70 

0.50 

1.18275E+00 

5.60423E+00 

3.29000E-01 

3.61259E-01 

9.8 

15.25 

0.30 

1.05379E+00 

4.27950E+00 

4.13000E-01 

3.75800E-01 

9.0 

25.40 

0.50 

1.18275E+00 

8.73365E+00 

2.71000E-01 

2.67972E-01 

1.1 
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a 


2a 

w 


INHERENT  FLAW  CONCEPT 
B/AL  (02/±45)s  EXPONENT  =  0.500 
UNNOTCHED  STRENGTH  =  800.1 


On 

Oo 


exp-1 


Co  = 

1.933202 

[ooj 

(ooJ 

TEST 

CALC 

%ERR 


0.25 

0.03 

1.00319E+00 

2.88986E-01 

8.78000E-01 

9.3801 6E-01 

6.8 

0.65 

0.07 

1.00787E+00 

6.0981 4E-01 

7.82000E-01 

8.58330E-01 

9.8 

1.25 

0.05 

1.00579E+00 

9.01570E-01 

7.21000E-01 

7.7481 5E-01 

7.5 

2.55 

0.05 

1.00590E-t-00 

1.94803E-t-00 

5.79000E-01 

6.5281 3E-01 

12.7 

2.55 

0.10 

1.011 51 E+00 

1.99769E+00 

5.71000E-01 

6.491 94E-01 

13.7 

5.10 

0.10 

1. 02251 E+00 

3.53966E-t-00 

4.64000E-01 

5.18313E-01 

11.7 

7.60 

0.30 

1.05343E+00 

4.77558E+00 

3.95000E-01 

4.27478E-01 

8.2 

12.70 

0.50 

1.18275E+00 

8.521 67E+00 

2.74000E-01 

3.0731  OE-01 

12.2 

15.25 

0.30 

1.05379E+00 

7.761 24E+00 

3.20600E-01 

3.18297E-01 

0.7 

25.40 

0.50 

1.18275E+00 

1.26673E+01 

2.28700E-01 

2.24854E-01 

1.7 

INHERENT  FLAW  CONCEPT 
B/AL  (0)6T  EXPONENT  =  0.500 
UNNOTCHED  STRENGTH  =  1672  Co  =  2.0801 1 4 


a 


w 


Y 


TEST 


/On' 

(oo^ 

CALC 


%ERR 


0.25 

0.03 

1.00319E+00 

4.85014E-01 

8.18000E-01 

9.41832E-01 

15.1 

0.65 

0.07 

1.00787E+00 

1.10108E+00 

6.64500E-01 

8.66060E-01 

26.5 

1.25 

0.05 

1.00579E+00 

1.82059E+00 

5.92000E-01 

7.85788E-01 

32.7 

2.55 

0.05 

1.00590E+00 

1.93787E+00 

5.80000E-01 

6.66334E-01 

14.9 

2.55 

0.10 

1.01 151  E+00 

2.17303E+00 

5.55000E-01 

6.62640E-01 

19.4 

5.10 

0.10 

1.01 151  E+00 

3.46239E+00 

4.68000E-01 

5.321 18E-01 

13.7 

7.60 

0.30 

1.05343E+00 

4.62366E+00 

4.00300E-01 

4.40046E-01 

9.9 

12.70 

0.50 

1.18275E+00 

6.44820E+00 

3.09800E-01 

3.17184E-01 

2.4 

15.25 

0.30 

1.05379E+00 

8.84757E+00 

3.02400E-01 

3.28768E-01 

8.7 

25.40 

0.50 

1.18275E+00 

1.22812E+01 

2.32000E-01 

2.3261 7E-01 

0.3 

C-2 


DISTRIBOTION  LIST 


No.  Of  Copies 


Kainan  Aerospace .  1 

Attn:  Mr.  A.  Falocxie 
Old  Windscxr  Road 
Blocmfleld,  CT  06002 

NAVAIWKVCEN . 1 

Tcm  Hess,  Code  60C1 

NAVAIRDE7CEN .  3 

Code  6064 

NAVAIRDEVCEN .  2 

Code  8131 

NAVAIRESVCEN . 25 

Code  6043 


DISmZBOTIGN  T.TfiT 


Mo.  Of  Copies 

Moclcwell  Intematicsial/Nbith  i^merican  Aircraft  Division .  1 

Attn:  Mr.  W.  O'Brien 

P.O.  BCQC  92098 

Los  Angeles,  CA  90009 

Rockwell  Internaticxial/MOrth  American  Aircraft  Division .  1 

Advanced  Technology  Program  Develcpnent 

Attn:  Dr.  ladkman 

P.O.  Box  92098 

Los  Angeles,  CA  90009 

Rockwell  InternatiOTal/Nbrth  American  Aircraft  Operatica^ .  1 

Attn:  Technical  Library 
P.O.  Box  582808 
Tulsa,  OK  74158 

Rockwell  InternatioTal/North  American  Aircraft  Operations .  1 

Attn:  Mr.  G.  Sherrick,  Engineering  Research  and  Technology 
P.O.  Box  582808 
TUlsa,  OK  74158 

Rockwell  Internaticaial/North  American  Aircraft  Operations .  1 

Attn:  Mr.  M.  Whitehead 
P.O.  Box  582808 
TUlsa,  OK  74158 

Rockwell  InternatiOTal/North  American  Aircraft  Operations .  1 

Attn:  Mr.  F.  Kaufinan 
P.O.  Box  582808 
TUlsa,  OK  74158 

Sikorslqr  Aircraft .  1 

Attn:  Technical  Library 
Nbrth  Main  Street 
Stratford,  CT  06601-1381 

Sikorsky  Aircraft .  1 

Attn:  Mr.  Samuel  P.  Garbo,  MS  S314A2 
North  Main  Street 
Stratford,  CT  06601-1381 

SiJcorsky  Aircraft .  1 

Attn:  Mr.  T.  Cook 
Nbrth  Main  Street 
Stratford,  CT  06601-1381 


DISTRIBOTIGM  LIST 


No.  Of  Copies 

McDonnell-Douglas  Astronautics .  1 

Attn:  Mr.  C.  M.  Roe 
5301  Bolsa  Avenue 
Huntingtcxi  Beach,  CA  92647 

Northrop  Aircraft  Corporation .  1 

Attn:  Technical  Library 
One  Northrop  Avenue 
Hawthorne,  CA  90250 

Northrop  Aircraft  Corporation .  1 

Attn:  Mr.  D.  Bold! 

One  Northrop  Avenue 
Hawthorne,  CA  90250 

Northrop  Aircraft  Corporation .  1 

Attn:  Dr.  M.  Ratwani 
One  Northrop  Avenue 
Hawthorne,  CA  90250 

Northrop  Aircraft  Corporation .  1 

Attn:  Mr.  B.  Butler 
One  Northrop  Avenue 
Hawthorne,  CA  90250 

Northrop  Aircraft  Corporation .  1 

Attn:  Mr.  R.  Whitehead,  Structures  Research 
One  Northrop  Avenue 
Hawthorne,  CA  90250 

Rockwell  Intematicavil/North  American  Aircraft  Division .  1 

Attn:  Technical  Library 

P.O.  Box  92098 

Los  Angeles,  CA  90009 

Rockwell  Intematicaial/North  American  Aircraft  Division .  1 

Attn:  Mr.  G.  Stewart,  Advanced  Structures  Engineering 

P.O.  Box  92098 

Los  Angeles,  CA  90009 

Rockwell  Internet ional/NOrth  American  Aircraft  Division .  1 

Attn:  Mr.  A.  Sanders 

P.O.  Box  92098 

Los  Angeles,  CA  90009 


DlfiTRIBUnOM  LIST 


McDormell-Dcuglas  Oorporatlon . . 

Douglas  Aircraft  Ocnpar^ 

Attn:  Mr.  A.  Hawley,  D^yt.  EIH,  IVS  212-10 
3855  Lakewood  Blvd. 
long  Beach,  CA  90846 

McDcnnell-Douglas  C3orporation . 

Douglas  Aircraft  Ocnpany 

Attn:  Mr.  R.  Petty,  Dept.  1L7,  IVS  35-98 

3855  lakewood  Blvd. 

long  Beach,  CA  90846 

McDonnell-Douglas  Corporation . . 

Douglas  Aircraft  Company 

Attn:  Mr.  Clifford  Y.  Kam,  Dept.  EL8,  212-20 
long  Beach,  CA  90846 

McDonnell-Douglas  Helicopter  Company . 

Attn:  Technical  Library 
Culver  City,  CA  90230 

McDonnell-Douglas  Heliccpter  Company . 

Attn:  Mr.  R.  M.  Verette 
culver  City,  CA  90230 

McDonnell-Douglas  Helicopter  Ccmpany . 

Attn:  Technical  Library 
5000  E.  McDowell  Road 
Mesa,  AZ  85205 

McDonnell-Douglas  Helicopter  Company . 

Attn:  Mr.  J.  K.  Sen 
338,  Bldg.  530 
5000  E.  McDowell  Road 
Mesa,  AZ  85205 

McDonnell-Douglas  Helicopter  Ccmpany . . . 

Attn:  Mr.  S.  Guymon,  MS  B337 
5000  E.  McDowell  Road 
Mesa,  AZ  85205 

Mc^)onnell-Douglas  Astronautics . 

Attn:  Technical  Library 
5301  Bolsa  Avenue 
HuntingtcMi  Beach,  CA  92647 


DISTRIBOnOM  LIST 


IHV  Aero^}ac3e  and  Defense  Ocnpany . . 

Vbui^t  Missiles  and  Advanced  Programs  Division 
Attn:  C.  Foreman,  Engineering  Project  Manager 
P.O.  Box  655907 
Dallas,  TX  75265-5907 

laV  Aero^^ace  and  Defense  Ocnpary . 

Vou^t  Missiles  and  Advanced  Programs  Division 
Attn:  Mr.  R.  Khi^t,  MS-ra-83 
P.O.  Box  650003 
Dallas,  TX  75265-0003 

lav  Aircraft  Product  Groc^ . 

Structures  Research  and  Development 
Attn:  Mr.  R.  Ely,  Technical  Project  Manager 
P.O.  Box  655907 
Dallas,  TX  75265-5907 

McDonnell-Douglas  Corporation . 

Attn:  Technical  Library 

P.O.  Box  516 

St.  Louis,  MD  63166 

McDonnell-Douglas  Corporation . 

Attn:  Mr.  C.  Pingle,  D^)t.  337,  IVS  0644324 

P.O.  Box  516 

St.  Louis,  MD  63166 

McDormell-Douglas  Corporation . 

Attn.  Mr.  D.  J.  Dorr,  D^.  390,  K/S  1021310 

P.O.  Box  516 

St.  Louis,  MD  63166 

McDonnell-Douglas  Oorporaticxi . 

Attn:  Mr.  C.  Saff,  Dept.  337,  H/S  0341160 

P.O.  Box  516 

St.  Louis,  MD  63166 

McDonnell-Douglas  Corporation . 

Attn:  Mr.  T.  Hinkel,  D^.  334,  IVS  0643063 

P.O.  Box  516 

St.  Louis,  MD  63166 

McDcxmell-Douglas  Corporation . 

Douglas  Aircraft  Ccnpany 
Attn:  Technical  Library 
Mail  Code  212-10 
long  Beach,  CA  90846 


DlgTRIBUTIOM  LIST 


lociUieed  Aeronaxitic  Systems  Ocnpany . 

Attn:  Mr.  E.  K.  Walker 

P.O.  Box  551  -  D^.  76-02,  B63,  Pit.  A1 

Burbank,  CA  91520 

Lockheed  Aeronautic  Systems  Oonpany . 

Attn:  C.  F.  Griffin,  D^jt.  76-24 
Advanced  Structures  Technology  D^jt. 
Burbank,  C3^  91520 

Loddieed-Georgia  Ccnpany . 

Attn:  Techniccil  Library 
86  South  Ccfcb  Drive 
Marietta,  GA  30063 

Loddieed-Georgia  Ocnpany . 

Attn:  Mr.  L.  R.  Meade,  MS  399 
86  South  Oobb  Drive 
Marietta,  GA  30063 

Lockheed-Georgia  Ocnpany . 

Attn:  Technical  Information 
D^)t.  72-34,  Zone  26 
Marietta,  GA  30063 

Lockheed-Missile  and  Space  Ocnpany . 

Attn:  Techniccil  Library 
1111  LocUieed  way 
Sunnyvale,  CA  94088 

Lockheed-Missile  and  Space  Ocnpany . 

Attn:  Mr.  J.  A.  Bailie 
P.O.  Box  504,  D^.  81-12,  B154 
1111  LocOdieed  Way 
Sunnyvale,  CA  94088 

mv  Aero^>aoe  and  Defense  Ocnpany . 

Vaught  Aercproducts  Division 
Attn:  Technical  Library 
P.O.  Box  655907 
Dallas,  TX  75265-5907 

mv  Aerospace  and  Defense  Ocnpany . 

Vou^t  Aercproducts  Divisicm 
Attn:  Mr.  J.  Piimi,  MS  194-51 
P.O.  Box  655907 
Dallas,  TX  75265-5907 


DISTRlBUnOM  LIST 


Mo.  of 


General  Electric  C3cii{)any; . 

Attn:  Dr.  B.  Eodini 
P.O.  Box  8555-M4018 
fhiladelphia,  PA  19101 

Oomiaan  Aeroqpaoe  corporation . 

Attn:  Technical  Library 
South  Oyster  Bay  Poad 
Bethpage,  Long  Island,  MY  11714 

Gruninan  Aero^>aoe  Corporation . 

Attn:  Mr.  R.  Hadoock,  Advanced  Materials  . 

South  Oyster  Bay  Road 

Bethpage,  Long  Island,  MY  11714 

Grumnan  Aerospace  Corporation . . . 

Attn:  Mr.  J.  Bruno,  B44-35 
South  Oyster  Bay  Road 
Bett^iage,  Ix»ig  Island,  MY  11714 

Grumman  Aero^»oe  Corporation . 

Attn:  Mr.  S.  Dastin 

South  Oyster  Bay  Road 

Betlpage,  Long  Island,  MY  11714 

Lockheed  Aeronautical  Systems  Ccxipany . 

Attn:  Technical  Library 
P.O.  Box  551 
Burbank,  CA  91520 

Lockheed  Aeronautical  Systems  Ocnpai^ . 

Attn:  C.  W.  Schneider 
D/70-03,  B^369,  P/B-6 
P.O.  Box  551 
Burbank,  CA  91520 

Locddieed  Aeronautic  Systems  Conpany . 

Attn:  Mr.  D.  E.  Pettit 
canyon  Research  Lab 
Burbank,  CA  91520 

Loddieed  Aeronautic  Systems  Ocnpany . 

Attn:  Mr.  R.  C.  Young 
Canyon  Research  Lab 
Buzbarik,  CA  91520 


OISISIBUTION  LIST 


No.  of 


General  Dynamics/Fort  Worth  Divisioi . 

Attn:  Mr.  T.  Kearns,  MS  6217 

P.O.  Box  748 

Fort  Worth,  TX  76101 

General  Dynamics/Fort  Worth  Dlvislcxi . . 

Attn:  Mr.  J.  A.  Fant 
P.O.  Box  748  -  MZ  2884 
Port  worth,  TX  76101 

General  Dynamics/Fort  Worth  Division . 

Attn:  Dr.  G.  Law,  Ocnposite  Structures  Certification 

P.O.  Box  748 

Fort  Worth,  TX  76101 

Genercil  Electric  Company . 

Attn:  Technic2d  Library 
1  Neumann  way 
Cincinnati,  <XI  45215 

General  Electric  Conpany . 

Attn:  Mr,  S.  Johnscai 
1  N^miann  Way 
Cincinnati,  OH  45215 

General  Electric  Ccmpany . 

Attn:  Mr.  A.  W.  Kim 
1  Neunonn  Way 
Cincinnati,  CH  45215 

General  Electric  Ccnpary . 

Attn:  Technical  Library 
P.O.  Box  8555 
Fhiladelphia,  PA  19101 

General  Electric  Ccnpary . 

Attn:  Mr.  A.  Garber 
i^.O.  Box  8555-M4018 
Fhiladeljhia,  PA  19101 

General  Electric  Ccnpary . 

Attn:  Mr.  C.  Zweben 
P.O.  Box  8555-M4018 
Hiiladelphia,  PA  19101 


DISTRZBUnON  LIST 


No.  of  Copies 

Boeing  Military  Aircraft  Gcnpany .  1 

Attn:  Mr.  J.  Avery 
P.O.  Box  20746 
WicJiita,  KS  67277-7730 

Boeing  Military  Aircraft  Ocnpaiiy .  1 

Attn:  Mr.  R.  Weu’ier 
P.O.  Box  20746 
Wichita,  KS  67277-7730 

Boeing  Military  Airplane  Ocnpany . . .  1 

Attn:  J.  luding,  MS  K3216 

P.O.  Box  20746 
Wichita,  KS  67277-7730 

Boeing  Military  Airplane  Ocnpany .  1 

Attn:  Mr.  R.  M.  Mayle,  MS  K06-72 
P.O.  Box  20746 
Wichita,  KS  67277-7730 

General  Dynamics/Oonvair  Division .  1 

Space  SystesQS  Divisim 
Attn:  Technical  LiJbrary 
P.O.  Box  85990 
San  Diego,  CA  92138 

General  Dynamics/Oonvair  Divisicai . 1 

Space  Systems  Divisicxi 
Attn:  Oraig  Rix,  MZ-23-8443 
P.O.  Box  85990 
San  Diego,  OA  92138 

General  Dynamics/Ocxwalr  Divisicxi .  1 

Space  Systems  Division 
Attn:  Mr.  D.  R.  Dunbar 
P.O.  Box  85990 
San  Diego,  OA  92138 


General  Dynamics/Fort  Worth  Divisicxi .  1 

Attn:  Technical  Library 

P.O.  Box  748 

Fort  Worth,  TX  76101 

General  Dynamics/Fort  Worth  Divisicxi .  1 

Attn:  Mr.  B.  Smith 

P.O.  Box  748 

Fort  Worth,  TX  76101 


DISTRIBOnCN  LIST 


NO.  of  copies 


Boeing  Helicxpter  Ocnpar^ .  1 

Attn:  Mr.  D.  Hart 
P.O.  Box  16858 
Diiladelphia,  PA  19142 

Boeing  Helicopter  Ocnpai^ .  1 

Attn:  Mr.  R.  L.  Pindkney 
P.O.  Box  16858 
Itdladelphia,  PA  19142 

Boeing  Helicopter  Oonpany .  1 

Attn:  Mr.  C.  Albrecht,  MS  P32-38 
P.O.  Box  16858 
fhiladelphia,  PA  19142 

^joeing  Helicopter  Oonpany . 1 

Attn:  Mr.  Brian  Lake,  MS  P30-18 
P.O.  Box  16858 
Ihiladelphia,  PA  19142 

Boeing  Oonpany .  1 

Attn:  Technical  Library 
P.O.  Box  3707 
SeatUe,  WA  98124-2207 

Boeing  Oonpany .  1 


Attn:  Mr.  R.  Hortcxi 
P.O.  Box  3707-je  3304 
Seattle,  WA  98124-2207 


Boeing  Ocnpary .  1 

Attn:  Mr.  M.  Oohodas 
P.O.  Box  3707,  MS  4A-16 
Seattle,  WA  98124-2207 

Boeing  Oonpany .  1 


Attn:  Mr.  J.  Quinlivan 
P.O.  Box  3707,  MS  6N-21 
Seattle,  WA  98124-2207 

Boeing  Military  Aircxaft  Oonpai^ 
Attn:  Technical  Library 
P.O.  Box  20746 
Wichita,  KS  67277-7730 


1 


DISIKIBUTICXI  LIST 


Mo.  of  Copies 

AVOO  Specialty  Materials  -  Textron . . .  1 

Attn:  Mr.  William  F.  (^^nt 
^peci2Q.ty  Materuds  Divisicn 
2  Industrial  Avenue 
louell,  MA  01851 

AVOO  Specialty  Materieds  -  Textron .  1 

Attn:  J.  Heroshaw 
Specialty  Matericds  Divisicxi 
2  Industrial  Avenue 
Lowell,  MA  01851 

Beech  Aircraft  Oorporaticxi .  1 

Attn:  Mr.  M.  B.  Goetz 
daypool  Bldg.,  Office  300 
4130  Linden  Avenue 
Dayton,  OH  45432 


Beech  Aircraft  Corporation .  1 

Attn:  Mr.  M.  P.  Djuric 
Kellogg  St.  and  W^h  Hoad 
Wichita,  KS  67201 

Bell  Helioopter/Textrcxi  Inc .  1 

Attn:  Technical  Lihrary 

P.O.  Box  482 

Port  Worth,  TX  76101 

Bell  Heliocpter/Textrcn  Inc .  1 

Attn:  Mr.  D.  Reisdorfer 

P.O.  Box  482 

Port  Worth,  TX  76101 

Bell  Helioopter/Textrcm  Inc .  1 

Attn:  Mr.  M.  K.  Stevenson 

P.O.  Box  482 

Port  Worth,  TX  76101 

Boeing  Helicopter  Oonpany .  1 


Boeing  Defense  and  ^)aoe  Groip 
Attn:  Dr.  C.  K.  Gunther 
P.O.  Box  16858,  MS  P30-30 
Philadelphia,  PA  19142-0858 

Boeing  Helicxpter  Oonpany . 

Attn:  Technical  library 
P.O.  Box  16858 
Philadelphia,  PA  19142 


1 


DISXRIBUnGM  LIST 


Mo.  Of 


C3cnnander . 

Navcil  Aviation  Depot 

Attn:  J.  Gresham,  Code  35210 

Mcas,  Cherry  Poiitt,  NC  28533-5030 

Ocninander. . . . . 

Naval  Aviation  Depot 

Attn:  M.  Lim,  Code  340 

NAS,  Jacksonville,  FL  32212-0016 

Ocmnander . 

Naval  Aviation  D^pot 

Attn:  J.  Deans,  Code  31310,  Bldg.  V88 

NAS,  Norfolk,  VA  23511-5899 

Ocninander . . 

Naval  Aviation  Depot 

Attn:  H.  Jarrett,  C3ode  32220 

NAS,  Norfolk,  VA  23511-5899 

Oommander . 

Naved  Aviaticxi  D^)ot 

Attn:  H.  Somraerfleck,  Code  36030 

NAS,  Norfolk,  VA  23511-5899 

Ocninanding  Officer . 

Naval  Aviation  Depot 
Attn:  D.  Perl,  Oode  343 

NAS,  North  Island,  San  Diego,  CA  92135-5112 

Oomnanding  Officer . 

Naval  Aviaticxi  D^xTt 

Attn:  M.  Williams,  Oode  31327 

NAS,  North  Island,  San  Diego,  CA  92135-5112 

Ocninanding  Officer . 

Naval  Aviation  Depot 
Attn:  D.  Yiapp,  Oode  342 
NAS,  Pensacola,  FL  32508-5300 

Ocmnander . . . 

Ikiited  States  Naved  Academy 
Attn:  Mechanical  Engineering  D^Tt. 
Annapolis,  MD  21402 


DlgntlBOTIOM  LIST 


No.  Of  Copies 

Ooninanding  Officer .  1 

U.S.  Army  R&T  laboratory  (ARRADOCXf) 

Attn:  K.  Abelsoi 
Building  182 
Dover,  NT  07801 

Ooonanding  Officer .  1 

U.S.  Any  R&T  laboratory  (ARRADOCN) 

Attn:  R.  Bank 
Building  182 
Dover,  NT  07801 

Ocninandiiig  Officer .  1 

U.S.  Any  ^plied  Technology  lab 
Attn:  G.  McAllister 
USARTL,  (AVRADQCXl) 

Fort  Eustis,  VA  23604-5418 

Ocnnending  Officer .  1 

U.S.  Any  ^jplied  Technology  lab 
Attn:  J.  Waller 
USARTL,  (AVRADOGM) 

Fort  Eustis,  VA  23604-5418 

NASA  Headquarters .  1 

Attn:  G.  Seidel 
CAST-Oode  RM 
washingtcxi,  D.C.  20546 


Administrator .  1 

National  Aeronautics  and  Space  Administration 
Attn:  Airframes  Branch,  FS  120 
Washington,  D.C.  20546 


Administrator .  1 

Rationed  Aeronautics  and  Space  Administration 
Langley  Research  Center 
Attn:  Dr.  C.  P.  Blankenship,  IVS-188M 
Hanpton,  VA  23365-5225 

Administrator .  l 

Rationed  Aeronenitics  and  ^)aoe  Administration 
langl^  Research  Center 
Attn:  Mr.  C.  E.  Harris,  MS  188E 
Hampton,  VA  23365-5225 


DISTSIBUnCXi  LIST 


No.  Of  Ocspies 


Oomnander .  1 

David  Taylor  Research  Center 
Attn:  A.  Macander,  Code  1720 
Bethesda,  MD  20084 

Ooranander .  1 

David  Taylor  Research  Center 


Attn:  E.  T.  Caitponeschi,  Code  2844 
Annc^lis,  MD  21402 


Ccratander .  1 

David  Taylor  Research  Center 
Attn:  R.  Crane,  Code  2844 
Annepolis,  MD  21402 

Ccmnander .  1 


David  Taylor  Research  Center 
Attn:  H.  Bdelstein,  Code  2870 
Arme^lis,  MD  21401 

Ccmnander . 1 

Naval  Surface  Weapons  Center 
Attn:  Dr.  J.  M.  Augl 
10901  New  Hanpshire  Avenue 
Silver  Spring,  MD  20903-5000 

Ccmnander .  1 

Naval  Surface  Veapars  Center 
Attn:  Dr.  J.  Goff,  R34 
White  Oak  Laboratory 
Silver  ^ring,  MD  20910 

Ccmnander . 1 

Naved  Aviation  Depot 
Attn:  N.  Amdur,  Code  342/N 
NAS,  Alameda,  CA  94501 


Ccmnander .  1 

Naval  Aviation  D^x>t 
Attn:  J.  Meyers,  Code  35110 
MCAS,  Cherry  Point,  NC  28533-5030 

Ccmnander . 

Naval  Aviation  Depot 
Attn:  J.  Riss,  Code  35430 
MCAS,  Cherry  Point,  NC  28533-5030 


1 


dishubutzon  list 


No.  of  Copies 


Ocmnamder .  1 

Naval  Air  Systems  Oomnand 
Attn;  AIR-53022 
Washington,  D.C.  20361 

Ocnnander .  1 


Naval  Air  Systems  Ocninand 
Attn;  AIRr53023 
Washington,  D.C.  20361 

Gcnmander . 1 

Naval  Air  Systems  Oommand 
Attn;  AIR-5304 
Washington,  D.C.  20361 

CQRinander .  1 

Naval  Air  Systems  Oonmand 
Attn;  AIR-53ar 
Washingtc»i,  D.C.  20361 

commander .  1 

Naval  Sea  Systems  command 
Attn;  Sea-05R 
Washington,  D.C.  20362-5101 

Gcnmander .  1 

Naval  Wea^xsns  Center 
Attn;  K.  Bailey,  Code  338 
China  lake,  CA  93355 


Conmander .  1 

Naval  Weapons  Center 

Attn;  J.  E.  Diamond,  Code  338 

China  LaJce,  CA  93355 

Conmanding  Officer . 1 

Naval  Ship  Engineering  Center 
Attn;  NSEC  6101E 
Arlington,  VA  20360 

Conmander .  1 

David  Taylor  Research  Center 
Attn;  R.  RodkNell,  Code  1720-6 
Bethesda,  ND  20084 


DISTRIBDnCN  LIST 


Nb.  of  Oopies 

Office  of  Naval  Technology .  1 

At±n:  W.  King,  afr-212 
800  N.  Quincy  Street 
Arlington,  VA  22217 

Office  of  Naval  Research .  1 

Attn:  A.  K.  Vasudevan,  Code  1216 
800  N.  Quincy  Street 
Arlington,  VA  22217 

Office  of  Navcd  Research .  1 

Attn:  Y.  Raj2^>a]cse,  Code  1132SM 
800  N.  Quincy  Street 
ArlingtcMi,  VA  22217 

Director .  1 

Naval  Research  Laboratory 
Attn:  Dr.  R.  Badciliance 
4555  Overlook  Avenue,  S.W. 

Washington,  D.C.  20375 

Director .  1 

Naval  Research  laboratory 
Attn:  Dr.  I.  Wolock,  Oode  6383 
4555  Overlook  Avenue,  S.W. 

Washington,  D.C.  20375-5000 

Director .  1 

Naval  Research  Laboratory 
Attn:  C.  Poranski,  Oode  6122 
4555  Overlook  Avenue,  S.W. 

Washington,  D.C.  20375-5000 

commander .  1 

Naval  Air  Systems  Ocnmand 
Attn:  AIR-530 
Washington,  D.C.  20361 


Ocmnander .  1 

NavcQ.  Air  Systems  OGranand 
Attn:  AIR-5302 
Washington,  D.C.  20361 


Ocmnander . 

Naval  Air  Systoans  Ocmnand 
Attn:  AIR-53021 
Washington,  D.C.  20361 


1 


DISTRIBOnON  LIST 


No.  Of  Ocpies 

Administxator . 1 

Naticxial  Aeronautics  and  Space  Administration 
George  C.  Marshall  Space  Fli^t  Center 
Attn:  Techniccd  Library 
Huntsville,  AL  35812 

Administrator .  1 

National  Aeronaxitics  and  Space  Administration 

Lewis  Research  Center 

Attn:  Dr.  C.  Chamis,  MS-49-6 

21000  Brooikpark  Road 

Cleveland,  OH  44153 

Administrator .  1 

National  Aeronautics  and  Space  Administration 

Lewis  Research  Center 

Attn:  M.  Hershberg,  MS  49-6 

21000  Broolqiark  Road 

Cleveland,  OH  44153 

Administrator .  1 

NatiCTial  Aeronautics  and  ^ce  Administration 

Lewis  Research  Center 

Attn:  Technical  Library 

21000  Broblqiark  Road 

Cleveland,  44153 

Administrator .  2 

Defense  Techniced.  Information  Center 
Bldg.  #5,  Oaineron  Station 
Alexandria,  VA  22314 


Federal  Aviation  Administration .  1 

Attn:  Mr.  J.  R.  Soderquist,  AIR-103 
800  Independence  Avenue,  S.W. 

Washington,  D.C.  20591 

Federal  Aviation  Administration. . . .  1 

Attn:  Mr.  L.  Neri,  Code  ACD-210 
Technical  Center 
Atlantic  City,  NT  08405 

Federal  Adviaticxi  Administration .  1 

Attn:  C.  Caiafa,  Code  ACr-330 
Technical  Center 
Atlantic  City,  MI  08405 


DISTREBOTIGM  LICT 


Mo.  Of 


Administrator . 

Naticsial  Aeronautics  and  Space  Administraticai 
Langley  Research  Center 
Attn;  Dr.  J.  Starnes,  MS  190 
Hanpton,  YA  23365-5225 

Administrator . 

Naticxial  Aerc»iautics  and  Space  Administraticai 
langl^  Research  Center 
Attn:  Dr.  M.  Mikulas,  Code  MS  190 
Hamptcxi,  VA  23365-5225 

Administrator . 

National  Aercxiautics  and  Space  Administration 
Langley  Research  Center 
Attn:  Mr.  W.  T.  Freeman,  MS  243 
Hampton,  VA  23365-5225 

Administrator . 

National  Aexx^iautics  and  Space  Administraticai 
Langley  Research  Center 
Attn:  Mr.  J.  W.  Deaton,  MS  188A 
Hamptcai,  VA  23365-5225 

Administrator . 

Nationcil  Aeronautics  and  Space  Administration 
langley  Research  Center 
Attn;  Mr.  J.  Davis,  MS  243/STPO 
Hampton,  VA  23365-5225 

Administrator . 

National  Aeronautics  and  Space  Administration 
l£mgley  Research  Center 
Attn:  Dr.  R.  Price 
Hanpton,  VA  23365-5225 

Adnunistrator . . 

Nationed  Aert^iautics  and  Space  Administration 
Langley  Research  Center 
Attn:  Dr.  G.  L.  Roderick 
Hampton,  VA  23365-5225’ 

Administrator . 

Naticmial  Aeronautics  and  Space  Administration 
George  C.  Marshall  Space  Fli^t  Center 
Attn:  R.  Schwin^iamer,  S&E-ASTN-M 
Huntsville,  AL  35812 


DISTRIBlfFIGM  LIST 


Mo.  of 


Ocninanding  Officer . 

Warner  Robbins  Ai"  logistics  OcamBnd 
Attn:  T.  F.  CJirLitian,  MBRD 
Robbins  Air  Force  Base,  GA  30198 

Ocninanding  Officer . 

Warner  Robbins  Air  logistics  Oonmand 
Attn:  W.  Schweinberg,  MMIRC 
Robbins  Air  Force  Base,  GA  30198 

Ocninanding  Officer . 

Picatinry  Arsenal 
PIASTEC 

Attn:  Librarian,  Oode  ERDRR-SCM-0,  Bldg.  351N 
Dover,  NJ  07801 

Ocninanding  Officer . 

U.S.  Army  Air  Mobility  R&D  Lab 
Attn:  H.  ReddicGc 
Fort  Eustis,  VA  23604 

Ocninanding  Officer . 

U.S.  Anty  Aviatic^  ^jplied  Technology  Directorate 

Attn:  T.  E.  Oondon 

SAVRT/IY-ASR 

Fort  Eustis,  VA  23604-5577 

Ocninanding  Officer . . 

U.S.  Amy  Aviation  ^plied  Technology  Directorate 

Attn:  Drew  Orlino 

SAVRT/TY-ATS 

Port  Eustis,  VA  23604-5577 

Ocninanding  Officer . 

U.S.  Army  Materials  and  Mechanical  Research  Oenter 
Attn:  D.  Oplinger,  SLCMT-MS 
Watertown,  MA  02171-0001 

Ocninanding  Officer . 

U.S.  Amy  Research  Office 
Durham,  NO  27701 

Ocninanding  Officer . 

U.S.  Anty  R&T  Laboratory  (AVRADOCM) 

Attn:  F.  Ininen,  Oode  DAVDL-AS-MS-207-5 
Ames  Research  Oenter 
Moffett  Field,  CA  94035 


DISTRIBOTIGM  LIST 


NO.  of 


Ocninandlng  Officer . . 

Research  and  Development  Center 
Attn:  FIBC/C.  Rams^ 

Wri^t  Battersc»i  Air  Force  Base,  CM  45433-6553 

CZcmnanding  Officer . . 

Wri^t  Research  and  Developnent  Center 

Attn:  FIBEC/Dr.  G.  Sendeciqrj 

Wri(^t  Pattersc«i  Air  Force  Base,  CH  45433-6553 

Qxtinanding  Officer . 

Wri<^t  Research  and  Developnent  CTenter 
Attn:  FIEH^.  F.  Wblfe 

Wriest  Patterson  Air  Force  Base,  (M  45433-6553 

ernmanding  Officer . 

Wri^t  Research  and  Developnent  Center 
Attn:  FIBA/L.  Kelly 

Wri^t  Patterson  Air  Force  Base,  CM  45433-6553 

ermranding  Officer . 

Wri^t  Research  and  Developnent  Center 
Attn:  MIHVDr.  J.  Whitney 

Wri^t  Patterson  Air  Force  Base,  CM  45433-6553 

ernmanding  Officer . 

Wri^t  Research  and  Develc^xnent  center 
Attn:  MLSEys.  Fecheck 

Wri^t  Patterson  Air  Force  Base,  CM  45433-6553 

Oomnanding  Officer . 

Wri^t  Research  and  Develcpnent  Center 
Attn:  MIHVDr.  M.  Knight 

Wri^t  Patterson  Air  Force  Base,  CM  45433-6553 

ceninanding  Officer . 

Wri^t  Research  and  Developnent  Center 
Attn:  MIB/F.  Cherry 

Wri^t  Patters(»i  Air  Force  Base,  CM  45433-6553 

Qannanding  Officer . 

Wti^t  Research  and  Developnent  Center 
Attn:  KLSE/T.  Reinhart 

wri^t  Patterson  Air  Force  Base,  CM  45433-6553 

Cenmanding  Officer . 

wriest  Research  and  Developnent  Center 
Attn:  MDIN/R.  C.  Tcinashot 

Wri^t  PattersCTi  Air  Force  Base,  CM  45433-6553 


DISIRlBirnCN  LIST 


No.  of  Copies 


Ocnnicuider .  1 

U.S.  Naval  Postgraduate  School 
Attn:  Professor  R.  Belli 
Mexiterey,  CA  93943 

Oenmander . . .  1 

U.S.  Naval  Postgraduate  School 


Attn:  Professor  M.  H.  Bank,  Code  67BP 
Monterey,  CA  93943 


C3cnitander .  1 

U.S.  Naval  Postgraduate  School 
Attn:  Professor  K.  Challenger 
Mcnterey,  CA  93943 

Oemnander .  1 

U.S.  Naval  Postgraduate  School 
Attn:  Technical  Library 
Mcaiterey,  CA  93943 

Oenmander . 1 

U.S.  Naved  Postgraduate  School 


Attn:  Dr.  E.  Rcbert  Wood,  Code  67 
Monterey,  CA  93943 

D^artment  of  the  Air  Force .  1 

Attn:  Dr.  M.  S2dkind 
Building  410 
Boiling  Air  Force  Base 
washingtcai,  D.C.  20331 

Department  of  the  Air  Force .  1 

Attn:  Dr.  J.  Amos 
Building  410 
Boiling  Air  Force  Base 
washingtoi,  D.C.  20331 

Oenmanding  Officer .  1 

Wri^t  Beseeuxh  and  Development  Center 
Attn:  FIBA/W.  Goesch 

Wri^t  Pattersen  Air  Force  Base,  OH  45433-6553 

Oenmanding  Officer . 

Wtight  Research  and  Develcpment  Center 
Attn:  FIB/R.  Bader 

Wright  Kitterson  Air  Force  Base,  OH  45433-6553 


1 


